TMB: Automatic Differentiation and Laplace Approximation

TMB is an open source R package that enables quick implementation of complex nonlinear random effect (latent variable) models in a manner similar to the established AD Model Builder package (ADMB, admb-project.org). In addition, it offers easy access to parallel computations. The user defines the joint likelihood for the data and the random effects as a C++ template function, while all the other operations are done in R; e.g., reading in the data. The package evaluates and maximizes the Laplace approximation of the marginal likelihood where the random effects are automatically integrated out. This approximation, and its derivatives, are obtained using automatic differentiation (up to order three) of the joint likelihood. The computations are designed to be fast for problems with many random effects (~10^6) and parameters (~10^3). Computation times using ADMB and TMB are compared on a suite of examples ranging from simple models to large spatial models where the random effects are a Gaussian random field. Speedups ranging from 1.5 to about 100 are obtained with increasing gains for large problems. The package and examples are available at http://tmb-project.org

The dynamics of social activation and suspiciousness in individuals at ultra-high risk for psychosis

INTRODUCTION: Social functioning is often impaired during the ultra-high risk (UHR) phase for psychosis, but group-level studies regarding the role of social functioning in transition to psychosis are inconsistent. Exploring the inter-individual differences which underlie the association between social functioning and psychotic symptoms in this phase could yield new insights.OBJECTIVE: To examine the idiographic and dynamic association between social activation and suspiciousness in individuals at UHR for psychosis using time-series analysis.METHODS: Twenty individuals at UHR for psychosis completed a diary application every evening for 90 days. Two items on social activation (quantity: 'time spent alone' and quality: 'feeling supported') and two items on suspiciousness ('feeling suspicious' and 'feeling disliked') were used. Time series (T = 90) of each individual were analyzed using vector auto regression analysis (VAR), to estimate the lagged (over 1 day) effect of social activation on suspiciousness, and vice versa, as well as their contemporaneous associations.RESULTS: Heterogeneous person-specific associations between social activation and suspiciousness were found in terms of strength, direction and temporal aspects.CONCLUSIONS: The association between social activation and suspiciousness differs amongst individuals who are at UHR for psychosis. These findings underline the importance of tailoring psychosocial interventions to the individual. Future studies may examine whether using results of single-subject studies in clinical practice to personalize treatment goals leads to better treatment outcomes.</p

Occupational exposure to gases/fumes and mineral dust affect DNA methylation levels of genes regulating expression

Many workers are daily exposed to occupational agents like gases/fumes, mineral dust or biological dust, which could induce adverse health effects. Epigenetic mechanisms, such as DNA methylation, have been suggested to play a role. We therefore aimed to identify differentially methylated regions (DMRs) upon occupational exposures in never-smokers and investigated if these DMRs associated with gene expression levels. To determine the effects of occupational exposures independent of smoking, 903 never-smokers of the LifeLines cohort study were included. We performed three genome-wide methylation analyses (Illumina 450 K), one per occupational exposure being gases/fumes, mineral dust and biological dust, using robust linear regression adjusted for appropriate confounders. DMRs were identified using comb-p in Python. Results were validated in the Rotterdam Study (233 never-smokers) and methylation-expression associations were assessed using Biobank-based Integrative Omics Study data (n = 2802). Of the total 21 significant DMRs, 14 DMRs were associated with gases/fumes and 7 with mineral dust. Three of these DMRs were associated with both exposures (RPLP1 and LINC02169 (2x)) and 11 DMRs were located within transcript start sites of gene expression regulating genes. We replicated two DMRs with gases/fumes (VTRNA2-1 and GNAS) and one with mineral dust (CCDC144NL). In addition, nine gases/fumes DMRs and six mineral dust DMRs significantly associated with gene expression levels. Our data suggest that occupational exposures may induce differential methylation of gene expression regulating genes and thereby may induce adverse health effects. Given the millions of workers that are exposed daily to occupational exposures, further studies on this epigenetic mechanism and health outcomes are warranted

Genetic variation in HSD17B13 reduces the risk of developing cirrhosis and hepatocellular carcinoma in alcohol misusers

BACKGROUND & AIMS: Carriage of rs738409:G in patatin-like phospholipase domain-containing 3 (PNPLA3) is associated with an increased risk for developing alcohol-related cirrhosis and hepatocellular carcinoma (HCC). Recently, rs72613567:TA in hydroxysteroid 17-beta dehydrogenase 13 (HSD17B13) was shown to be associated with a reduced risk for developing alcohol-related liver disease and to attenuate the risk associated with PNPLA3 rs738409:G. This study explores the risk-associations between these two genetic variants and the development of alcohol-related cirrhosis and HCC. APPROACH AND RESULTS: Variants in HSD17B13 and PNPLA3 were genotyped in 6,171 participants, including: 1,031 with alcohol-related cirrhosis and HCC; 1,653 with alcohol-related cirrhosis without HCC; 2,588 alcohol misusers with no liver disease; and 899 healthy controls. Genetic associations with the risks for alcohol-related cirrhosis and HCC were determined using logistic regression analysis. Carriage of HSD17B13 rs72613567:TA was associated with a lower risk for both cirrhosis (OR 0.79 [95% CI 0.72-0.88], p=8.13×10-6) and HCC (OR 0.77 [95% CI 0.68-0.89], p=2.27×10-4), while carriage of PNPLA3 rs738409:G was associated with an increased risk for developing cirrhosis (OR 1.70 [95% CI 1.54-1.88], p=1.52x10-26) and HCC (OR 1.77 [95% CI 1.58-1.98], p=2.31×10-23). These associations remained significant after adjusting for age, sex, body mass index, type II diabetes mellitus and country. Carriage of HSD17B13 rs72613567:TA attenuated the risk for developing cirrhosis associated with PNPLA3 rs738409:G in both men and women but the protective effect against the subsequent development of HCC was only observed in men (p=1.72×10-4; ORallelic, 0.75; 95% CI, 0.64-0.87). CONCLUSIONS: Carriage of variants in PNPLA3 and HSD17B13 differentially affect the risk for developing advanced alcohol-related liver disease. A genotypic/phenotypic risk score might facilitate earlier diagnosis of HCC in this population

Second intravenous immunoglobulin dose in patients with Guillain-Barre syndrome with poor prognosis (SID-GBS):a double-blind, randomised, placebo-controlled trial

Background Treatment with one standard dose (2 g/kg) of intravenous immunoglobulin is insufficient in a proportion of patients with severe Guillain-Barre syndrome. Worldwide, around 25% of patients severely affected with the syndrome are given a second intravenous immunoglobulin dose (SID), although it has not been proven effective. We aimed to investigate whether a SID is effective in patients with Guillain-Barre syndrome with a predicted poor outcome. Methods In this randomised, double-blind, placebo-controlled trial (SID-GBS), we included patients (>= 12 years) with Guillain-Barre syndrome admitted to one of 59 participating hospitals in the Netherlands. Patients were included on the first day of standard intravenous immunoglobulin treatment (2 g/kg over 5 days). Only patients with a poor prognosis (score of >= 6) according to the modified Erasmus Guillain-Barre syndrome Outcome Score were randomly assigned, via block randomisation stratified by centre, to SID (2 g/kg over 5 days) or to placebo, 7-9 days after inclusion. Patients, outcome adjudicators, monitors, and the steering committee were masked to treatment allocation. The primary outcome measure was the Guillain-Barre syndrome disability score 4 weeks after inclusion. All patients in whom allocated trial medication was started were included in the modified intention-to-treat analysis. Findings Between Feb 16, 2010, and June 5, 2018, 327 of 339 patients assessed for eligibility were included. 112 had a poor prognosis. Of those, 93 patients with a poor prognosis were included in the modified intention-to-treat analysis: 49 (53%) received SID and 44 (47%) received placebo. The adjusted common odds ratio for improvement on the Guillain-Barre syndrome disability score at 4 weeks was 1.4 (95% CI 0.6-3.3; p=0.45). Patients given SID had more serious adverse events (35% vs 16% in the first 30 days), including thromboembolic events, than those in the placebo group. Four patients died in the intervention group (13-24 weeks after randomisation). Interpretation Our study does not provide evidence that patients with Guillain-Barre syndrome with a poor prognosis benefit from a second intravenous immunoglobulin course; moreover, it entails a risk of serious adverse events. Therefore, a second intravenous immunoglobulin course should not be considered for treatment of Guillain-Barre syndrome because of a poor prognosis. The results indicate the need for treatment trials with other immune modulators in patients severely affected by Guillain-Barre syndrome. Funding Prinses Beatrix Spierfonds and Sanquin Plasma Products. Copyright (C) 2021 Elsevier Ltd. All rights reserved

Controlling bias and inflation in epigenome- and transcriptome-wide association studies using the empirical null distribution

We show that epigenome- and transcriptome-wide association studies (EWAS and TWAS) are prone to significant inflation and bias of test statistics, an unrecognized phenomenon introducing spurious findings if left unaddressed. Neither GWAS-based methodology nor state-of-the-art confounder adjustment methods completely remove bias and inflation. We propose a Bayesian method to control bias and inflation in EWAS and TWAS based on estimation of the empirical null distribution. Using simulations and real data, we demonstrate that our method maximizes power while properly controlling the false positive rate. We illustrate the utility of our method in large-scale EWAS and TWAS meta-analyses of age and smoking

Dutch juvenile idiopathic arthritis patients, carers and clinicians create a research agenda together following the James Lind Alliance method: A study protocol

Background: Research on Juvenile Idiopathic Arthritis (JIA) should support patients, caregivers/parents (carers) and clinicians to make important decisions in the consulting room and eventually to improve the lives of patients with JIA. Thus far these end-users of JIA-research have rarely been involved in the prioritisation of future research. Main body: Dutch organisations of patients, carers and clinicians will collaboratively develop a research agenda for JIA, following the James Lind Alliance (JLA) methodology. In a 'Priority Setting Partnership' (PSP), they will gradually establish a top 10 list of the most important unanswered research questions for JIA. In this process the input from clinicians, patients and their carers will be equally valued. Additionally, focus groups will be organised to involve young people with JIA. The involvement of all contributors will be monitored and evaluated. In this manner, the project will contribute to the growing body of literature on how to involve young people in agenda setting in a meaningful way. Conclusion: A JIA research agenda established through the JLA method and thus co-created by patients, carers and clinicians will inform researchers and research funders about the most important research questions for JIA. This will lead to research that really matters

Mendelian randomization integrating GWAS and eQTL data reveals genetic determinants of complex and clinical traits

Genome-wide association studies (GWAS) have identified thousands of variants associated with complex traits, but their biological interpretation often remains unclear. Most of these variants overlap with expression QTLs, indicating their potential involvement in regulation of gene expression. Here, we propose a transcriptome-wide summary statistics-based Mendelian Randomization approach (TWMR) that uses multiple SNPs as instruments and multiple gene expression traits as exposures, simultaneously. Applied to 43 human phenotypes, it uncovers 3,913 putatively causal gene-trait associations, 36% of which have no genome-wide significant SNP nearby in previous GWAS. Using independent association summary statistics, we find that the majority of these loci were missed by GWAS due to power issues. Noteworthy among these links is educational attainment-associated BSCL2, known to carry mutations leading to a Mendelian form of encephalopathy. We also find pleiotropic causal effects suggestive of mechanistic connections. TWMR better accounts for pleiotropy and has the potential to identify biological mechanisms underlying complex traits